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Amateur-Radio Applications of the 


Fast Fourier Transform 


Part 1 


This article is more than just a description of 
the latest software for the DSP Computer, as 
interested readers will soon find out. Never- 
theless, it is also a continuation of the series 
of articles about DSP Techniques and written 
by the same author. 


1, 
INTRODUCTION 


The Fast Fourier Transform (FFT) is a very 
efficient numerical algorithm to compute a 
discrete Fourier transform. The FFT algorithm 
has found many different applications and it 
is not limited only to digital computers and digital 
signal processing. For example, a completely 
analog implementation of the FFT algorithm can 
be made with a number of 3 dB directional 
couplers and delay lines to feed an antenna array. 


Obviously the FFT algorithm is widely used in 
digital signal processing too. In this article a 
particular OSP application of the FFT algorithm 


is described and discussed in detail: FFT spec- 
trum analysis. Since the FFT algorithm can be 
implemented on almost any digital computer and 
a FFT spectrum analyzer only requires a little 
additional hardware, FFT spectrum analysis is 
very convenient for many practical applications. 


In order to understand the operation of a FFT 
spectrum analyzer and its advantages and draw- 
backs when compared with other spectrum 
analysis techniques, a short description of the 
various Fourier transforms and the operation of 
the FFT algorithm is included. Building a practi- 
cally working spectrum analyzer around the bare 
FFT algorithm is also discussed. No mathe- 
matical proofs are given just to keep the discus- 
sion as simple as possible. The interested reader 
can find the former in almost any book describing 
digital signal processing. 


Although FFT spectrum analysis is presently 
limited to the audio frequency range, at least for 
amateur resources, it has many amateur-radio 
applications ranging from weak signal detection 
to modulation analysis. A few typical amateur 
applications are discussed later in this article, 
including spectrum function plots and intensity 
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spectrograms obtained from real! signals using 
a FFT spectrum analyzer. 

Finally, a FFT spectrum-analyzer program for 
the DSP computer, described in UKW-BE- 
RICHTE/VHF COMMUNICATIONS, is pre- 


saith tassiah ‘ok Soak: Malay ales ak he 
computer screen on a laser printer. 


2. 

THE FOURIER TRANSFORM AND 
THE DISCRETE FOURIER 
TRANSFORM 


The Fourier transform is a mathematical opera- 
tion that computes a new function F(w) from an 
origina! function {(t), as shown in fig. 2.1. Both 
functions have real arguments t and w while their 
values f(t) and F(w) are in general complex 
numbers. F(w) is called the Fourier transform of a 


trom a given F(w) is very similar to the Fourier 


transform itself. 
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The Fourier transform is frequently used in 
physics to compute the frequency spectrum of a 


pressure, 
electrical current or voltage) as a function of time. 
In all physical problems f(t) is a real function. The 
new variable w represents the frequency and the 
function F(w) is the frequency spectrum of the 
physical quantity observed. The frequency 
spectrum is a complex function of a real variable 
«@, The absolute value of F(w) represents the 
magnitude of a given spectral component and the 
argument of F(w) represents the phase relative 
to the (chosen) time origin. Since f(t) is always 
& real function one does not really need to com- 
pute F(w) for negative frequencies: F(—w) is 
simply the complex-conjugate of F(w) for a real 
A(t). 


While the Fourier transform is a powerful analyti- 
cal tool to perform theoretical computations 
it has at least two constraints which could never 
be fulfilled in practice: infinite bandwidth and 
infinite frequency resolution (which implies an 
infinite observation time), regardies of the 
method, analog or digital, used to perform the 
Fourier transform. 


at least twice the signal bandwidth. The Fourier 


F(w) 


Fig. 2.1.: 
Fourier Transform 
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transform of a sampled signal is shown in fig. 2.2.. 
The integral is replaced with a sum, the bounds of 
the sum are, however, still infinite. An important 
difference should be noticed in plot of F(w): the 
spectrum of a sampled signal is a periodic 
function and its period is inversely proportional 
to the sampling period. It is therefore sufficient 
to compute (either in an analog or in a digital 
way) just one period of F(w), usually between 
~ nldtand + n/At. 


The (original) Fourier transform integral (or sum) 
has infinite bounds: the integration (or summa- 
tion) should be performed from minus infinity to 
plus infinity. Of course no real-world signal will 
ever last that long! It is therefore completely 
sufficient to compute the integral or the sum only 
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over the time interval when the signal exists. inte- 
grating over a finite amount of time limits the 
frequency resolution, which is inversely pro- 
portional to the integration time. This implies 
that the frequency spectrum could also be re- 
presented by discrete samples in place of a 
continuous function. 

The procedure that computes a finite number 
of spectral lines from a finite number of signal 
samples is called the Discrete Fourier Transform 
(DFT) and is shown in fig. 2.3. To properly 
describe the signal spectrum one needs at least 
the same number of frequency samples as there 
are input signal samples. Of course, to obtain 
meaningful results, the frequency interval (spec- 
tral line spacing) has to be chosen in a close 
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relationship to the time interval (signal sampling 
step). Similar constraints also apply to a com- 
pletely analog Fourier transtorm (analog spec- 
trum analyzer). 


In order to simplify computations, both time and 
frequency units are usually normalized to 1 as 
shown In fig. 2.4. Time t now only takes the 
values of 0, 1, 2,... (N — 1) and frequency w also 
takes the values of 0, 1, 2,...(N — 1), ifthe discrete 
Fourier transform produces N spectral lines from 
N signal samples. As a result of this normalisation 
the constant 2 «7/N appears in the complex ex- 
ponent function. This constant is chosen such 
that the resulting spectral lines cover exactly one 
period of the periodic signa! spectrum. 
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Fig. 2.4.: 
Normalizing time and 
frequency units 
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The DFT can be. computed on any general- 
purpose computer. Its operation is similar to N 
bandpass FIR filters, each filter having N stages 
and tuned to its own frequency. The DFT is not 
computationally efficient, since the number of 
computations required increases with N’. If the 
OFT is computed on a real signal (real f(t)) with N 
signal samples, then the result only includes 
N/2 spectral lines ranging from zero to half the 
sampling frequency. The remaining N/2 spectral 
lines are simply complex conjugates that can be 
computed in a much simpler way once the first 
N/2 spectral lines are known, reducing the total 
number of computations required to about N’/2. 


Will be continued 
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Amateur-Radio Applications 
of the Fast Fourier Transform 


Part2a 


3. 
THE FAST FOURIER TRANSFORM 
(FFT) ALGORITHM 


The Discrete Fourier Transform supplies a very 
useful result, bul unfortunately requires a very 
large number of computations on a digital com- 
puter or a very large number of components if 
performed by an analog circuit. A more efficient 
algorithm that provides exactly the same result as 
DFT with a considerably smaller effon required is 
called the Fast Fourier Transform. As an ex- 
ample, to compute a 1024 data point DFT the FFT 
algorithm requires about 200 times less compu- 
tations (or analog components) than a straight- 
forward DFT. Funher, the number of compu- 
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lations required in the FFT algorithm is only 
proportional to N/2 « log 2 (N). The FFT algorithm 
can therefore be performed on a very large 
number of data samples without significantly 
increasing the number of computations per 
sample. 


The basic building block of the FFT algorithm is 
called a “butterfly operation. A “butterily” 
operation consists of a phase-shift operation and 
a sum/difference operation. It operates on two 
input variables and produces two results. A single 
“putterfly“ operation can aiready compute a 2- 
point FFT, as shown in fig. 3.1. In the case of a 
2-point FFT, the two input variables are the input 
data !o the FFT algorithm ano the two results are 
already the result of the algorithm. The phase 
shift is equal to 7 in this case. 


p F(1) = (O0)e" + ((1)0” 


Fig. 3.1.: 
Two-point FFT — a single 
“butterfly* operation 
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Fig. 3.2. shows how a 4-point FFT works. A 4- 
point FFT is computed in two stages. Each stage 
includes two "butterfly" operations. In the first 
Stage, all phase shifts are equal to 7. In the 
second stage, the phase shifts are a/2 and z. 
Note that each output of the first stage “butter- 
flies” is fed to exactly one input of the second 
stage “butlerflies’. Considering the periodicity of 
the complex exponent function, the four ouipuls 
correspond exactly to the result obtained with a 
siraightiorward DFT, the latter, however, requires 
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Fig. 3.2.: 

Four-polnt FFT, computed 
each using two “butterlly“ 
operations 
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16 phase shifts (8 if one does nat consider zero 
phase shifis) compared to the 4 phase shifts of 
the FFT algorithm. 


Similarly, a 8-point FFT 1s computed in three 
stages as shown in fig. 3.3. Each stage includes 
4 “butterlly” operations. Again, in the first stage 
all phase shifts are equal tor. In (he second stage 
ihe phase shifts are 7/2 and 7. \n the third stage 
the phase shifts are 7/4, 7/2, 3 « 1/4 and r. The 
algorithm block diagram follows a regular pattern 
too, suggesting that a FFT algorithm working on 
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Fig. 3.3.: Elght-point FFT, computed in three steps, each using four “butterfly” operations 
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an arbitrary large number of samples could be 
designed in a similar way. In the case of an 8-point 
FFT, only 12 “butterflies” are required compared 
to the 64 shif/add operations of a straightlonward 
DFT algorithm to get the same result. 


To design FFT algoritms, operating on an even 
larger number of data samples, one should 
therefore investigate the passibility of combining 
several FFTs computed on a smaller number of 
samples. Fig. 3.4. shows how to combine two 
N-point transforms (not necessarily computed 
using the FFT algorithm) into one single transform 
on 2 N (twice the number of) points. In addition 
to the two N-point transforms, N “butierllies“ 
are required. These “butterflies” require N dif- 
ferent phase shifts ranging from n/N, 2 = a/N, 
3 « a/N...(N — 1) * 7/N, 7 in steps of 7/N. 


The principle shown in fig. 3.4. is in fact used to 
Gesign a FFT algorithm operating on any data 
length N that is a power of 2. A FFT algorithm can 
thus operate on 2, 4,8, 16, 32. 64, 128, 256, 512, 
1024... dala points. The 4point FFT can be 
derived from the 2-poim FFT, the 8-point FFT 
can be derivec from the 4-point FFT, the 16-point 
FFT can be derived from the 8-point FFT etc... 
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N * butterflies” Fig. 3.4.: 
o F(a) Obtaining a doubte- 
p F(2) length (2N) FFT from 
two single-length (N) 
FFTs 


Each doubling of the data poinis only requires 
an addninal stage so the total number of stages 
is equal to log 2 (N), Each stage requires N/2 
“butterlies*. The total number of “butterfly” 
operations is therefore equal to N/2 « log 2 (N). 


A quick Look at fig. 3.3. shows that the results 
do not appear in any reasonable order at the 
ouiput of the FFT algorithm: 1, 5, 3, 7, 2,6, 4, 0. 
Comsidenng the construction principle shown in 
fig. 3.4., a very simple rule to find the desired 
outpul can be found. This rule is called bit- 
reversed addressing and it is shown in fig. 3.5. 
for the 8-poim FFT example. To find an output 
number one has to take the corresponding input 
number, write this number down in binary format, 
reverse the order of bits, conven the number 
back to decimal and add 1. This rule can be 
easily implemented in digital hardware (dedi- 
cated OSP microprocessors), where dala is 
stored in consecutive memory locations. 


To implement the FFT algorithm, a number of 
“butterily” operations have to be periormed. 
An analog implementation called the “Butler 
mainx” uses delay lines to obtain phase shifts 
and “ral-race“ hybnds or 3 dB directional 
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couplers for the add/subiract operations. The 
inputs of the FFT circuit are connected to the 
elements of a linear antenna array and the out- 
puts of the FFT circuit to receivers and/or trans- 
mitters. Since the required phases and magni- 
tudes of the signals feeding an antenna array 
are the Fourier transform of the desired radiation 
pattern, the outputs of the FFT circuit correspond 
directly to beams in the various directions. 


Implementing the FFT algorithm on a digital 
computer most operations are performed with 
complex numbers. Although phase shifts are 
easier to perform if complex numbers are held 
in @ magnitude/phase Jormat, sums and differ- 
ences require a_ real/imaginary-component 
number format. Since conversions from one 
number format to another are very time-con- 
suming, all computations are usually done in 
the real/fimaginary-component format. In this 
case a phase-shift operation requires four 
teal multiplications with coefficients trom a 
precomputed table of phase shifis and two 
additions. Each summation of two complex 
numbers requires two rea} summations and 
each complex difference requires two real 
differences. 


The phase-shift coefficients are precomputed 
and stored in memory, since only N different 
coefficients are required in all stages of a N-point 
FFT. The same coefficients can be used in the 
following FFT, if a number of FFTs have to be 
computed on changing data. Finally, the phase- 
coefficient table Includes sines and cosines which 
are much more time-consuming (o compute than 
the multiplications and addnions required for the 
FFT itself. 


\f a FFT is computed on real data, then only half 
of the outputs contain interesting data, the other 
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half are just complex conjugates. To use the 
algorithm more efficiently, another set of input 
Gata can be fitied into the imaginary part of the 
input variables. After ihe FFT algorithm is per- 
formed, the two results can be separated by 
simple additions and subtractions using the 
symmetry laws of the Fourier transform. If de- 
sired, ihe two results can be further combined 
inte a single, double length FFT. 


The inverse DFT can easily be performed using 
the FFT algorithm in the reverse direction. The 
humber of mathematical operations is identical 
except for an additional division by N for each 
data point to obtain the original magnitude back. 


There are even more efficient algorithms to 
compute a DFT or its inverse. All of ihem are, 
however, based on the FFT principle described 
above and require a little more programming 
efforts to further reduce the number of com- 
pulautons required. 


4. 
SPECTRUM ANALYSIS USING THE 
FFT ALGORITHM 


One of tha most obvious applications of the FFT 
algorithm is a FFT spectrum analyzer. The FFT 
algorithm itself is performed on a digital computer, 
usually a DSP microprocessor. The input signal is 
Provided in a digital format {rom an A/D convener. 
The microprocessor itself can display the result in 
a variety of formats. 


However, to build a digital FFT spectrum analyzer 
some additional functions are required. A block 


£(0) >>---> 0 = OOOB >>---> 000B = 0 O + 1 = 1 >>---> F(1) 
£(1) >>---> 1 = 001B >»---> 100B = 4 4 +41 = 5 >>---> F(5) 
£(2) >»>---> 2 = 010B >>---> 010B = 2 2+ 1 = 3 >>---> F(3) 
£(3) >>---> 3 = 011B >>---> 110B = 6 6 + 1 = 7 >>---> F(7) 
£(4) >>---> 4 = 100B >>---> 001B = 1 1+ 21 = 2 >>---> F(2) 
£(5) >>---> 5S = 101B >>---> 101B = 5 S +12 6 >>---> F(6) 
£(6) >>---> 6 = 110B >>---> 0D11B = 3 3 + 1 = 4 >>---> F(4) 
£(7) >>---> 7 = 111B >>---> 111B = 7 7+ 1 = B >>---> F(O) 
Fig. 3.5.: Output bit-reversed mapping for a eight-polnt FFT 
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Fig. 4.1.: Block diagram of an FFT spectrum analyzer 


diagram of a DSP microprocessor-based spec- 
trum analyzer is shown in fig. 4.1. The analog 
input signal is first sent to an analog lowpass 
(bandpass) filter to prevent aliasing like in any 
DSP application. The analog filler is followed 
by a sample-and-hold and an A/D converter, The 
A/D converter feeds a buffer memory, since the 
FFT algorithm operates on blocks of data. 
Before the FFT the signal samples may be 
weighted optionally. The complex spectrum pro- 
vided by the FFT algorithm is used to compute the 
magnitudes of the spectral components, the 
phase information is discarded, Averaging is 
used ta improve the signal-to-noise ratio in 
some measurements. An optional linear-to- 
logarithmic conversion is standard for all spec- 
trum analyzers. Finally, the result can be dis- 
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display 


played in different formats and a conventional 
computer printer can be used to obtain a hard- 
copy. 


The performance of a DSP FFT spectrum ana- 
lyzer is mainly limited by the performance of the 
A/O converter used. The conversion speed of the 
A/D converter defines the maximum bandwidth 
and the resolution of the A/D converter defines 
the available dynamic range. After the A/D 
converter the cigitized input signal can be con- 
veniently stored in mermary if the microprocessor 
1s unable to process the data in real time. Also 
all suitable microprocessors offer a computa- 
tional accuracy of at least 16 - 24 bits allowing a 
much wider dynamic range than any A/D con- 
verter, 
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Fig. 4.2,: 
Overlapping FFTs 
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Fig. 4.3.: Weighting signal samples 


Even if the DSP microprocessor is fast enough 
to process all of the dala in real time, the input 
data stream has to be sent to a buffer memory 
first, since the FFT algorithm operates on blocks 
of data, not on single samples. Without weighting 
FFTs are usually performed on contiguous, but 
non-overlapping blocks of data to obtain almost 
all of the spectral information contained in the 
input signal. 


In the case of input signal weighting, the con- 
inbution of the samples at the beginning or at the 
end of a block is very limited. To use all of the 
information contained in the input signa) over- 
lapped FFTs have to be performed. Both cases 
ara shown in fig. 4.2. In the second case the 
overlap is sel to about 50 % for a raised-cosine 
weighting function. 


\f the FFT algorithm is performed on a raw, 
unweighted block of data, the corresponding 
spectrum is distorted. Since the FFT does not 
consider any data samples before the block nor 
any data samples after the processed block of 
data, Ihe actual input signal to the FFT algorithm 
corresponds to the real input signal modulated 
with a rectangular impulse of the length of the 
data block. A rectangular impulse with steep 


DERLEING “TS times Pu 
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SVEPLCING 


Fig. 4.4.: No weighting (above) versus raised- 
cosine welght\ng (below), same input 
signal, LOG vertical scale 


leading anc trailing edges has a very wide (re- 
quency spectrum of the form sin(X)/X. The 
resulting output will be the real signal spectrum 
convoluted with the sin(X)/X function. 


Although the real signal spectrum can not be 
abiained since it requires an infinite amount 
af time, a much more accurate spectrum can be 
obtained by weighting the signal samples as 
shown in fig. 4.3. Weighting means multiplying 
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each signa! sample with a constant whose value 
depends on the position of the sample inside the 
input data block. The practical effect of weighting 
is to replace the abrupt ON/OFF transitions with 
smooth transitions at the beginning and at the and 
of the data block. Weighting functions are se- 
lected !o minimize the distortion of the signal 
spectrum. Raised-cosine and Gaussian functions 
are usual choices, since they have a very narrow 
own spectrum with low side lobes. 


FUNCTION FLOT = AVERAGING d tires 


Fig. 4.5.: No averaging (above), averaging 4 times 
(below), same input signal and other 
settings ag in fig. 4.6. 
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The efiects of no weighting versus weighting 
are shown in fig. 4.4. The two plots show the 
spectrum of the same signal obtained in two 
different ways. A linear 512-point frequency 
scale is used on the horitzontal axis and a 
logarithmic amplinude scale (15 bits or 90 dB/ 
full scale) is used on the vertical axis. All the 
parameters, including the input signal, are the 
same for both plots except for weighting: the 
plot above was obtained without any weighting 


My 
FUNCTION PLOT: GVERAGING =~ L6 times Aud 
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Flg. 4.6.: Averaging 16 times (above) or 64 times 
(below), same input signal and other 
settings as in fig. 4.5. 
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while the plot below was oblaines with a raised- 
cosine weighting. Raised-cosine weighting 
clearly provides much more clear spectral lines, 
with no INEXISTENT sidebands. On the other 
hand, weighting slightly reduces the frequency 
resolution making the peaks broader. Selecting 
weighting or not is therefore a tradeoti between 
dynamic range and frequency resolution. 


lf a measurement is affected by random noise. 
averaging among a number of otherwise 1dentical 
measurements improves the accuracy of the 
resull. Automatic averaging is very easy to imple- 
ment on any microprocessor-controlled test 
equipment. Averaging is used in spectrum 
analyzers 10 improve the signal-to-noise ralio of 
the displayed data. The improvement that can 
be obtained by averaging is shown in fig. 4.5. 
and fig. 4.6. All four plois were obtained from the 
same input signal in an identical way except lor 
different amounts of averaging. 


The lineas-to-logarithmic conversion includes 
the computation of a Jogarihm for each spectral 
component. Since the logarithm ts a rather *slow* 
function on digital computers, a look-up-lable 
algornhm or a similar approach has to be used 
fo avoid unnecessary loading of ihe computer. 
The same constraint applies to the display 
procedure: most computers require more lime 
to draw a high-resolution plot than to compute 
the FFT algorithm. The display routine and/or 
dedicated hardware has to be quick enough to 
avoid slowing-down the spectrum analyzer. On 
the other hand, a hardcopy of the video display 
ig usually vary easy to obtain on any computer 
using a standard printer or plotter, al least when 
compared to analog instrumentation with CRT 
displays. 


Besides the conventional frequency/amplitude 
function plot other types of display are possible 
on digital computers. A practically very useful 
type of display is the intensity spectrogram. In 
the latter frequency is still plotted on the hori- 
zontal axis. Each FFT result is, however, re- 
presented by a single-image line and the pixel 
brightness is used to represent the magnitude 
of a speciral component. The results of succes- 
sive FFTs are plotted on successive lines, 
showing the results of a large number of consecu- 
tive measurements on jus! One Compuler screen. 
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Intensity spectrograms are useful when analyzing 
continuously changing signals 


The display of a FFT spectrum analyzer usually 
includes the complete frequency range covered 
by the FFT algorithm: from zero to half the signal 
sampling frequency The frequency resolution is 
then simply equal to the frequency span divided 
by the number of lines displayed. The frequency 
resolulion may be slightly worse if weighting is 
used in front of the FFT algorithm. Of course 
partial displays are possible too, showing just a 
subset of speciral lines computed by the FFT 
algorithm. The sample-and-hotd circuit in front 
of the A/D convener is an excellent harmonic 
mixer Ii is therefore possible to observe a dif- 
ferent frequency band just by replacing the input 
lowpass filter with a bandpass filler for the 
selected frequency range. 


Finally, a companson has to be made between 
a FFT-based spectrum analyzer and a scanning- 
receiver type spectrum analyzer (conventional 
analog RF spectrum analyzer) Of course a 
scanning-receiver type specirum analyzer could 
be implemented on a digital computer as well. 
A FFT-based spectrum analyzer has, however, 
a very imponant advantage over a scanning- 
receiver spectrum analyzer: regardless of the 
hardware used the FFT spectrum analyzer uses 
ihe available spectral information in a much 
more efficient way resulting in a much quicker 
operation. 


As an example, consider that a 5 kHz wide fre- 
quency band has to be analyzed to a resolulion 
of 10 Hz. A scanning receiver with a 10 Hz bano- 
width has to dwell on each 10 Hz frequency step 
for about 0 1 seconds, resulling in a total sweep 
time of about 50 seconds! On the other hang, a 
FFT-based spectrum analyzer needs to sample 
a 5 kHz wide signal wilh a sampling frequency 
of 10 kHz A 1024-point FFT has to be used to 
obtain a 512-point display, so the total "scanning" 
time is 0.1024 seconds! 


In the above real-world example the FFT spec- 
\rum analyzer is abou) 500 times faster! The 
reason for (his ts (hat a conventional scanning- 
receiver specirum analyzer only uses the in- 
formation contained in ils receiver bandwidth, all 
the other informabon contained in the signal 
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is simply rejected! On the other hand, the FFT- 
based spectrum analyzer uses all of the in- 
formation contained in the signal since the 
FFT algonihm corresponds to a bank of 512 
parallel bandpass filters in the above example. 
Such a bank of filters would be prohibitively 
expensive and difficult to make using conven- 
tional analog technolagy. 


A FFT-based spectrum analyzer can therefore 
be used in applications where a conventional 
scanning-receiver type spectrum analyzer 1s not 
practical due to the tao long scanning time or 
completely useless since the signal is not 
available for ihe scanning time period required. 
fven in the case when the scanning-receiver 
type can be used, the FFT-type can provide a 
much more accurate result in the same ume, 
averaging among a large number of measure- 
ments. 


Unfortunately the bandwidth and dynamic range 
of digital FFT spectrum analyzers are severely 
limited by the available hardware, mainly A/D 
converters. li is therefore necessary to under- 
siand the advantages and disadvantages both 
techniques to select ihe most suitable one for 
@ panicular problem, since the two techniques 
are complementing each other rather than 
competing al the present state of technology. 


Will be continued! 
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Fast Fourier Transform 


Part 2b (Concluding) 


5. 
AMATEUR APPLICATIONS OF AFFT 
SPECTRUM ANALYZER 


Although FFT spectrum analysis is presently 
limited to the audio-frequency range or slightly 
above, it has many interesting and very useful 
amateur-radio applications. In the following 
section a few typical amateur-radio applications 
will be presented including the spectrum plots 
and intensity spectrograms obtained. All of the 
latter were obtained by connecting the output of 
an amateur SSB or FM receiver to the MC68010 
based DSP computer described in a series of 
articles in UKW-BERICHTE/VHF COMMUNI- 
CATIONS. All of the plots were obtained with 


A FFT spectrum analyzer is a useful tool when 
building a SSB receiver or transceiver. One of 


the most difficult tasks when building a SSB 
receiver is to measure the passband of the 
crystal filters used, regardless whether are they 
homebrew or commercially available items. To 
obtain a reliable result, a very stable sweep 
generator is required in addition to a storage 
oscilloscope due to the slow sweeping speed 
required, Alternatively, a FFT spectrum analyzer 
can be connected to the receiver audio output 
and a wideband noise source to the receiver 
input (if the receiver own noise is not sufficient). 
Thanks to the speed of the FFT spectrum analysis 
the result can be obtained quicker than with the 
sweep generator, allowing real-time tuning of the 
trimmers in and around the crystal filters. Two 
typical results, plotted on a logarithmic amplitude 
scale, are shown on fig. 5.1.: a 2 KHz SSB filter 
above and a 500 Hz CW filter below. Using FFT 
spectrum analysis, the tradeoff between accu- 
racy and speed is selected by choosing the 
averaging factor. However, even with no aver- 
aging, the FFT analysis will only provide a noisy 
plot while a too fast sweep generator will provide 
a completely distorted and thus useless result. 

A FFT spectrum analyzer is able to reliably 
detect very weak signals hidden in noise, far 
beyond what a human ear can do, since it is not 
limited to certain frequency bands, resolution 


219 


VHF COMMUNICATIONS 4/90 


Fig, 5.1.: Audio frequency response of a SSB receiver: left with SSB filter, right with CW filter 


bandwidths or averaging intervals. The plots 
on fig. 5.2. were obtained by tuning a 2m amateur 
SSB receiver to a weak CW beacon and then 
decreasing the signal level with an input attenu- 
ator to obtain a signal-to-noise ratio of about 
—15 dB in a 2.5 kHz bandwidth or —5 dB ina 
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250 Hz bandwidth. Although the signal-to-noise 
ratio is around +10 dB in the spectrum analyzer 
resolution bandwidth, some averaging is required 
to reliably detect the signal. On the frequency/ 
amplitude plot (above on fig, 5.2.) the averaging 
was performed by the computer and then the 
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Fig. 5.2.: Weak signal detection; averaging, LOG scale (left) versus intensity spectrogram (right) 
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Fig. 5.3,; CW signal with a slight “FM-ing" (left) and affected by a quickly changing Doppler shift (right) 


result was plotted. On the intensity spectrogram 
(below on fig. 5.2.) no averaging was performed 
by the computer. Averaging is, however, per- 
formed by our eyes when observing the spectro- 
gram! 

Even weaker signals could be detected either 
by increasing the frequency resolution or by 
increasind the averaging factor or both. The 
practical limit is mainly imposed by the time 
required for the signal to be available for a reliable 
detection. In a practical communications system 
there are other constraints too: receiver and 
transmitter frequency instability or phase noise 
and propagation effects. The FFT spectrum 
analyzer can solve the problem of frequency 
uncertainty, since It allows to observe a frequency 
band of a few kHz instantaneously with a reso- 
lution of 5 to 10 Hz. The phase noise of trans- 
mitters and receivers should be minimized any- 
way. Unfortunately, some propagation effects 
also show up as phase or amplitude noise, 
especially in EME (moonbounce) communica- 
tions. These effects are proportional to the carrier 
frequency so major advantages of using FFT 
techniques for EME communications can only be 
expected on VHF and UHF frequencies. On 
these frequency bands FFT signal processing 
may decrease the RF link performance re- 


quirements by 10 to 20 dB, since higher figures 
would result in useless data rates. In any case, 
a communications protocol has to be agreed 
upon before these techniques can be used: hand- 
keyed CW is certainly not a good choice for 
computer processing, 


Fig. 5.3. shows how a Morse-keyed CW trans- 
mission looks on an intensity spectrogram. The 
keying was too quick for the dots and dashes to 
appear on the spectrogram, the interruptions 
correspond to the longer spaces between letters. 
In order to be able to see the single dots, much 
more frequent FFTs should be made and each 
FFT should be taken on less data samples since 
time and frequency resolution are of course 
reciprocal. 

Receiver AGC effects are easily visible, in- 
creasing the noise level during longer pauses. 
The spectrogram above was obtained from the 
AO-13 Mode-B 2 m beacon and shows a slight 
“FM-ing"; the trace has “tails” to the right at the 
beginning of each transmission. The spectrogram 
below was obtained from the LUSAT-1 70 com CW 
beacon. A quickly changing Doppler frequency 
shift is easily visible, as well as an interference 
to the left corresponding to the third harmonic 
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SPECTROGROM | RESOLUTION 2200 He / full scale STOP 


Fig. 5.4.: FSK RTTY signal spectrum: LOG vertical scale (left) and Intensity spectrogram (right) 


of the signal, generated somewhere in the audio 
stages of the SSB receiver. 


A FFT spectrum analyzer can be used to identify, 
tune-in and measure the parameters of a FSK 
RTTY transmission as shown on fig. 5.4. The 
two RTTY tones are easily visible both on the 
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amplitude/frequency display and on the intensity 
spectrogram. Both of them were obtained from 
the AO-13 50 baud RTTY Mode-B beacon. In 
the middie of the spectrogram a period with no 
keying is visible: only one of the tones, unmodu- 
lated, is transmitted during this period. During 
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Fig. 5.5.: AFSK 1200 bps packet-radio signal spectrum: LOG vertical scale (left) and spectrogram (right) 
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Fig. 5.6.: 1200 bps satellite PSK packet-radio transmission: flags (left) and data (right) 


active keying the tone traces become wider due 
to the 50 bps modulation sidebands. 


A FFT spectrum analyzer can therefore be con- 


Both tone frequencies and the keying shift can 
be quickly and accurately determined. in addition, 
a FFT display provides information about signal 
distortion (selective fading) or interferences. 
Especially in the latter case a FFT display will 
provide some useful information about the 
countermeasures to be taken (notch filter tuning 
for example) and about their effectiveness. 


The parameters of a FSK or AFSK signal are 
more difficult to identity if the data rate is com- 
parable to the frequency shift, as shown on 
fig. 5.5. 1200 bps packet-radio uses 1200 Hz 
and 2200 Hz tones, the shift is therefore 1000 Hz, 
comparable to the data rate of 1200 bps. The 
spectrum of such a transmission is an almost 
contiguous frequency band with just a few peaks, 
that do not necessarily correspond to the tone 
frequencies as shown on fig, 5.5, above. The 
intensity spectrogram below shows the inter- 
mittent nature of packet-radio signals. Due to the 


A FFT spectrum analyzer can be used as a 
valuable tuning aid for PSK packet-radio com- 
munications. Phase-shift keying is used for 


The spectrum of a 1200 bps packet-radio trans- 
mission from the PACSAT-1 satellite is shown on 
fig. 5.6. The plot above shows the spectrum 
during the transmission of flags between the 
packets. On this plot it is easy to identify the 
carrier, surrounded by sidebands spaced at 
150 Hz. During the transmission of packets con- 
taining random data the spectrum looks almost 
like perfect noise (fig. 5.6. below) and it is much 
more difficult to identify the correct tuning, In 
fact a straightforward PSK transmission contains 
very little redundancy and there are no residual 
carriers either. 


Some more redundancy can be noticed in the 
AO-13 400 bps BPSK Mode-B beacon trans- 
mission due to manchester coding. Fig. 5.7. 
above shows the spectrum during the trans- 
mission of the filling bytes (50 H) in between the 
data frames: this repetitive pattern generates 
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Fig. 5.7.: AO-13 400 bps BPSK telemetry spectrum: LOG vertical scale (left) and intensity spectrogram (right) 


the many discrete spectral lines. On the intensity 
spectrogram on fig. 5.7. below it is easy to identify 
the single data frames and the filling byte periods 
in-between them, Even data frames sometimes 
contain repetitive patterns causing discrete 
spectral lines inside the data frames. Regardiess 
of the modulation data the signal spectrum has 
two main lobes caused by the manchester 
coding. The latter are positioned symmetrically 
around the carrier, suppressed by this modulation 
technique. 


A FFT audio spectrum analyzer is useful to check 
WEFAX and APT satellite signals and related 
receiving equipment. Using a FFT spectrum 
analyzer important details of an unknown satellite 
transmission, like that of a new satellite, can 
easily be determined. Further, the receiving 
equipment can be checked for a correct de- 
emphasis (to avoid loosing geometrical reso- 
lution) and sources of eventual interferences. 


All this is easier to describe on a well known 
example, like the Meteosat WEF AX transmission 
shown on fig. 5.8. On the frequency/amplitude 
plot on fig. 5.8. above the strongest spectral 
component is the 2400 Hz subcarrier. The sub- 


224 


carrier should be surrounded by symmetrical 
sidebands If the deemphasis is properly adjusted. 
The sidebands depend of course on the picture 
content. The most notable detail are the two 
symmetrical peaks corresponding to the line- 
sync 840 Hz bursts. These have a rather rounded 
peak, since their duration is very short: their 
spectral width is inversely proportional to their 
duration. The intensity spectrogram shown on 
fig. 5.8. below shows that the sync bursts appear 
in about every second FFT conversion. Also, the 
spectrum close to the 2400 Hz subcarrier also 
changes with the picture line period. The anno- 
tation transmitted at the end of the picture is 
well visible, followed by the discrete spectral 
lines of the 450 Hz stop tone lasting 5 seconds. 
What follows is just part of the spectrum of a DCP 
retransmission between two WEFAX pictures. 
The DCP flags create two strong discrete spectral 
line traces on the spectrogram, interleaved with 
three DCP data frames visible on fig. 5.8. below, 


The spectrum of a NOAA APT ‘transmission, 
shown on fig. 5.9. is a little more complex. There 
still a strong 2400 Hz subcarrier component, but 
there are two different sync bursts of 832 Hz and 
1040 Hz respectively. The sidelobes generated 
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Fig, 5.8.; METEOSAT WEFAX signal spectrum: LOG vertical scale (left) and intensity spectrogram (right) 


by both sync bursts are well visible on fig. 5.9. Although a FFT spectrum analyzer is actually 
above. On the intensity spectrogram on fig. 5.9. an audio-frequency test equipment, it has many 
below it is easy to notice that these sync pulses _—interesting applications in the amateur-radio 
appear interleaved in the APT signal. The re- _ field. The above examples were chosen just to 
maining spectrum also changes with the same _— show a few possible applications of a FFT spec- 
period. trum analyzer and the many different ways a 
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Fig. 5.9.: NOAA-10 APT signal spectrum: LOG vertical scale (left) and intensity spectrogram (right) 
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FFT spectrum analyzer can be used, Probably 
most of the applications are yet to be discovered, 
since FFT spectrum analysis only became 
available to amateurs with recently developed 
low cost and high performance microprocessors 
and A/D converters. 


6. 

THE FFT SPECTRUM-ANALYZER 
PROGRAM FOR THE DSP 
COMPUTER 


The FFT algorithm can be implemented on almost 


stored in the same two memory locations where 
the arguments were taken from. Therefore the 
memory requirements are very low, including 
just a buffer of the size of the input data block 
and another buffer of a similar size to hold the 
phase coefficients. 


Of course the execution time will depend on the 
computer used. Computing a 1024-point complex 


FFT on an IBM AT compatible equipped with 
the mathematical coprocessor and executing a 
compiled program will still require several sec- 
onds. A well-written machine code program for a 
general-purpose 16-bit microprocessor can 
compute a 1024-point FFT in a few hundred 
milliseconds using 16-bit integer arithmetics. 
Finally, dedicated DSP microprocessors are 
able to compute a 1024-point FFT in tens of 
milliseconds with the top-of-line devices re- 
quiring a few milliseconds. 


AFFT spectrum-analyzer program was written for 
the DSP computer published in UKW-BE- 
RICHTE/VHF-COMMUNICATIONS. This com- 
puter includes a 86-bit logarithmic A/D converter 
input port and a 512-pixel by 256-line, 256-grey- 
level video display, The MC68010 microproces- 


sor used in this computer is able 
a 1024-point complex FFT in about 


The FFT spectrum-analyzer program for the DSP 


commands is shown on fig. 6.1. The same HELP 
menu is obtained each time a wrong command is 
lines for the function plot or spectrogram and 
the 8 bottom lines to annotate the main program 
parameters. 

The program parameters are selected by typing 
their initial letter as shown on the HELP menu. 
The selected parameter will then appear in the 
annotation text line. The selected parameter 
value can be then adjusted: increased or de- 
creased, by using respectively the + and — keys. 
The * key will set ALL of the parameters to their 
default values. 


The program parameters include; 

AVERAGING: averaging factor, ranging from 1 
(no averaging) to 4096 and adjustable in powers 
of 2. 


BORDER, DIVISIONS, ERASER, INK, PAPER: 
grey-level count, adjustable between 0 and 248 
in steps of 8. 
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program showing all implemented 
commands 
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HARDCOPY: each time H is depressed, a hard- 
copy file PLOT.DAT is created. The format is 
selected to be understood by most dot-matrix 
and laser printers. Since these printers can not 
print grey levels, the treshold between black 
and white is set to a grey-level count of 127/128! 


OVERLAP: switches the overlap (50 %) option 
ON or OFF (toggle operation). When ON, the 
overlap option is indicated by an “O" in the 
annotation text. 


RESOLUTION: selects the sampling frequency 
and therefore the full-scale frequency span or 
resolution between 301 and 6400 Hz by setting 
the sampling-frequency divider modulo, 


SPECTROGRAM DISPLAY; selects an intensity 
spectrogram type display, as indicated in the 
annotation text. 


TRIGGER; starts or stops (toggle operation) 
the operation of the spectrum analyzer. The 
status RUN or STOP is indicated in the annotation 
text. In the STOP mode, the display is “frozen”. 


Mt / ful 


Fig. 6.2.: TP 3040 input filter passband (left) and MK 5156 A/D converter noise with no input signal (right) 
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selects a LOG video scale (default), 


WEIGHTING: switches the raised-cosine 
weighting option ON or OFF (toggle operation). 
When ON, the weighting option is indicated by 
a “W" in the annotation text. 


The performance of this FFT spectrum analyzer is 
of course limited by the analog hardware used: 


simply connecting a wideband noise source to 
the input. The dynamic range of the 8-bit logarith- 
mic A/D converter used is only about 40 dB due 
to quantisation noise. The performance of the 
spectrum analyzer is however better: in many 
practical cases the display dynamic range 
reaches 60 dB since the quantisation noise 
usually appears as wideband noise. 


The input switched-capacitor filter and A/D con- 
verter own noise is weil visible on fig. 6.2. below. 
Using a logarithmic amplitude display the quanti- 
sation steps become very large at low signal 
levels. In the LOG vertical mode, the full scale 
amplitude range is 15 bits or 90 dB. The dif- 
ference between counts 1 and 2 should there- 
fore correspond to the difference between 
counts 16384 and 32768 on a LOG scale, how- 
ever in the former case there are no additional 
counts possible between 1 and 2 resulting in 
relatively large “steps” on the display, 

With a 10 MHz CPU clock, the FFT spectrum- 


The input data buffer routine is designed to skip 
data blocks automatically if the CPU Is unable 
to further process them in real time, like in the 
case of a higher sampling frequency or over- 
lapping FFTs. Sometimes it is very useful to be 
able to operate in these conditions even if some 
data is lost. On the other hand, an analog scan- 
ning-receiver type spectrum analyzer usually 
discards over 99 % of the information contained 
in the input signal! 


The FFT program includes a 1024-point complex 


The FFT spectrum-analyzer program is written 
partially in the MC68010 machine code and 


Finally, it is hoped that this article will encourage 
amateurs to use the FFT algorithm and FFT 
spectrum analysis, since the latter can easily be 
implemented on the described DSP computer 


reason for this is probably the poor understanding 
of this new technique by the majority of spectrum- 
analyzer users. 
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CX 120P 
Art.No.: 0500 DM 43.50 


suitable for PCB installation 

Freq. range 500 MHz; Pmax 150 W; insertion loss 
< 0.2 dB/1 GHz; isolation 40 dB/1 GHz; SWR 1.06/ 
1 GHz; 11 - 16 VDC; 120 mA/13.8 VOC 


Cx 140D 
. Art.No.: 0501 DM 59.50 


2x cable connect. RGSBU, 1 N-socket 

Freq. range 2.5 GHz (max); Prax 200 W/500 MHz; 
insertion loss 0.2 dB/2.5 GHz; Isolation 40 dB/ 
1 GHz; SWR 1.06/1 GHz; 11 - 16 VDC; 120 mA’ 
13,8 VOC 


CX 230 
Art.No.: 0502 DM 82.50 


SS 
3 x BNC sockets 
Freq. range 1.5 GHz; Pmax 300 W/500 MHz; inser- 
tion loss 0.2 dB/500 MHz; isolation 30 dB/500 MHz; 
SWR 1,1/500 MHz; 11 - 15 VDOC; 160 mA/12 VDOC 


CX 520D 
Art.No.: 0503 DM 99.50 
3x N-sockets, with free-contact grounding 
Freq. range 2.5 GHz; Pmax 300 W/1 GHz; inser- 
tion loss 0,2 dB/1.5 GHz; isolation 50 dB/1 GHz; 
SWRA 1.1/1 GHz; 11 - 15 VDC; 160 mA/12 VOC 


2 CX 600 N 
Art.No.: 0504 DM 89.— 
3x N-sockets 
Freq. range 1.5 GHz; Pmax 600 W/S00 MHz; in- 
sertion loss 0.2 dB/500 MHz; isolation 35 dB/500 
MHz; SWR 1,1/500 MHz; 10 - 15 VDC; 160 mA/ 
12 VDC 


CX 600 NC 
Art.No.: 0505 DM 84.— 
3x cable connect. RG 58, 1 N-socket 

Freq. range 1 GHz; Pmax 600 W/500 MHz; in- 
sertion loss 0.2 dB/500 MHz; Isolation 35 dB/500 
MHz; SWR 1.1/500 MHz; 10 - 15 VDOC; 160 mA/ 
12 VOC 
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